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Abstract 

Many fluctuating systems consist of macroscopic structures in addition to noisy signals. Thus, 
for this class of fluctuating systems, the scaling behaviors are very complicated. Such phenom- 
ena are quite commonly observed in Nature, ranging from physics, chemistry, geophysics, even to 
molecular biology and physiology. In this paper, we take an extensive analytical study on the "gen- 
eralized detrended fluctuation analysis" method. For continuous fluctuating systems in arbitrary 
dimensions, we not only derive the explicit and exact expression of macroscopic structures, but also 
obtain the exact relations between the detrended variance functions and the correlation function. 
Besides, we undertake a general scaling analysis, applicable for this class of fluctuating systems 
in any dimensions. Finally, as an application, we discuss some important examples in interfacial 
super roughening phenomena. 
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I. INTRODUCTION 



Fluctuating systems have been a study of great interest for their generic scaling behaviors 
widespread in Nature [ll Much attention has recently focused on interfacial superrough- 
ening phenomena for their anomalous scaling behaviors, different from the conventional 
scaling ansatz|3|]. Many experiments have been performed, including growth of nickel sur- 
faces by pulse-current electrodepositionjj], cultivated brain tumor growth in a Petri dishful, 
spontaneous imbibition of a viscous fluid by a model porous mediump, and molecular 
beam epitaxial growth j?], ^. All these superrough interfacial growth processes consist of 
local interfacial orientational instability. The boundary conditions (due to the finite size 
of real experiments) restrict the infinite development of orientational instability Thus, the 
interfacial configuration gradually develops global mountains or valleys as time increases. 
The appearance of macroscopic structures is the signature of superrough interfacial growth 
processes. Interestingly, in a different context of science (computational molecular biology) , 
with some specific mapping rules the DNA sequences can be viewed as a random walkP.llf|. 
Lots of data show that the DNA sequences form mosaic structures (patchiness) . Peng et 

nnn 

al. 9 : , HO, IUI| have proposed a "generalized detrended fluctuation analysis" method to de- 
termine the long-range correlation in DNA sequences but exclude the effect caused by the 
mosaic structure. The basic idea of this method is to eliminate the trend (or the macroscopic 
structure), which is expressed in terms of a polynomial with the coefficients determined by 
numerical least-squares- fit. It has been adopted as a standard analysis method in various 



time series problems, for example, the heart-rate 



s between neural action potentials 
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fluctuations during sleep jl^j], interfiring 



, and temporal fluctuations in seismic 



time interva ; 
sequencesjl4|. 

However, few analytical works have been done regarding the "generalized detrended fluc- 
tuation analysis" method. Since many fluctuating systems are modeled by stochastic partial 
differential equations, the explicit expression of the macroscopic structures (or the trends) 
is necessary for further analytical study of such systems. The numerical fitting can only 
give the values of various coefficients but not their relations with the statistical quantities 
of the systems. In addition, when a high-order detrending is necessary, the numerical least- 
squares-fit of many coefficients consumes large amounts of computation time. Besides, what 
are the explicit relations between the statistical quantities of the systems before detrending 
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and after detrending? Can a general scaling analysis be derived? Thus, we are strongly mo- 
tivated to take an analytical study on the fluctuating systems with formation of macroscopic 
structures. In Ref. jl5j . we have studied the fluctuating systems in one dimension. In this 
paper, we plan to work out all these issues for fluctuating systems in arbitrary dimensions. 
Moreover, we notice that many superrough interfacial growth processes consist of a very 
distinct property: the local roughness exponent equal to 1, regardless of other exponents or 
parameters of the systems. What is its origin? What does it imply? Are all these problems 
related? In this paper, we also plan to explore this intriguing phenomenon in some detail. 

The outline of this paper is as follows. In Sec. II, we present the core ideas of the "gen- 
eralized detrended fluctuation analysis" method. Then, we derive the exact expression of 
macroscopic structures in fluctuating systems of arbitrary dimensions, expressed in terms of 
the Legendre polynomials. In Sec. Ill, the exact relations between the correlation function 
and the detrended variance functions are derived. In Sec. IV, through the detailed investiga- 
tion of the kernel functions, we give a general scaling analysis applicable for any fluctuating 
systems with macroscopic structure formation. In Sec. V, as an application, we focus on 
some important examples in interfacial superroughening processes. In the last section, a 
brief summary is given. 



II. SYSTEMATIC EXTRACTION OF MACROSCOPIC STRUCTURES 



In this section, we will explicitly derive the exact expression of the macroscopic structure 
of fluctuating systems in arbitrary dimensions. First, let y(x) denote a random function of 
x, with x being a vector in the <i-dimensional space. One of the most important statistical 
quantities is the variance of ?/(x) in an observation window: 



w 



\i=l 



<v(*)> 



(i) 



with (. . ^-r-rd , _ denoting the average over an observation window, centered at x and of size 
Ilf=i h- The corresponding statistical quantity, related to macroscopic structure extraction, 
is the {g}-th ({q} = {q 1 , . . . , q d }) degree detrended variance of y(x) in that observation 
window: 



\i=l 
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where the coefficients {a{ n }} of the {g}-th degree macroscopic structure polynomial ?%}(x; x) 
[= Yl\n}={o} a {n} Hi=i( x i ~ %i) ni with the sum running over the set of indices < rti < for 
all 1 < i < d] are determined by dw 2 q y/da^ n } = 0. Hence, the coefficients {ct{ n }} satisfy the 
following relation 



M 



i=l '- 1 In'WO) i=l 



'h 



rii+n'- 



2(m + n' + 1) 



(3) 



{n'}={0} 

To determine {a{ n }} ; one needs to deal with the inverse of a tensor of rank 2d. By choosing 
a set of orthonormal functions as the basis, we can actually avoid the difficulty in solving 
the inverse of a high-rank tensor. 

Let {f n (x);n = 1,2,...} denote such a set of orthonormal functions on [—1,1] (i.e., 
(/m(^)/m'(^))[-i,i] = <Wn')- 2%}( x ;x) is then expressed as 



M d 
y w (x;x) = J2 c mI1/« 

{„}={0} i=l 



(4) 



The requirement <9w /<9c{ n } = leads to 



C{n} 



y(x) n A 



2 = 1 



^2[Xi Xi 



nti 



(5) 



Therefore, 



{?} / d 

y {?} (x;x) = ^ (2/(x)n/n 

{n}={0} \ 



8=1 



U 



(6) 



The above result can be viewed as the original y(x) convolved by some specific filter. Based 
on the problems to be tackled, one can choose a specific set of orthonormal functions to 
facilitate their analysis. In this paper, we plan to analyze the scaling behaviors of the 
correlation function and the detrended variances of the systems. Thus, the set of Legendre 
polynomials 



[i/2] 



^) = E(-i; 



• (2z-2j-l)!! 



-x 



{ - 2j for i = 1,2,..., 



is a perfect candidate to be chosen. From the orthogonal property 

Oii! 



(7) 



(Pi{x)P v {x)) 



2i + l 



(8) 
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we then define our basis by 

fi(x) = y/2i + lPi(x). (9) 

Besides, in the following analysis we will use the symbol [q] to denote the union of the 
sets of indices {n} = {ni, n 2 , ■ ■ ■ , n^} with |{n}| = n\ + n 2 + • • • + rid < Q and all rii being 
nonnegative integers. Hence, 



y[?](x;x)= £ 

IWI<9 \ 



i=l 



) n + 1)^ 



(10) 



As an illustration, let us consider the case with d = 2. With the dimensionless quantity 
3j = (xi — Xi)/li, we have the coefficients 



(y( X ))/lX/ 2 ;x 

2v / 3(y(x)2;i) ilXi2; x 



coo 

ClO 

Cll = 12(y(x)ziZ 2 )i lX Z 2 ;x 

C20 
C21 
C30 



(y(x)(l2^-l)> IlXla , 



15<3/(x)(12z? - l)z 2 ) IlX , a;5 
v / 7(y(x)(20^-3^ 1 )), lXi2;x 



and the corresponding y[0](x; x) to y[3](x; x) are also displayed as follows 
y[0](x;x) = (s/(x)), lX , a;i 

2 

y[l](x;x) = (y(x))j lX , a;x + 12^^(y(x)^ lX j a;x 

i=i 

2 r5 -i 
y[2](x;x) = y[l](x;x)+£ -(12z 2 - l)(y(x)(12z 2 - l)), lX i 2 ;x + 144^ 2 (y(x)^ 2 ) il>d2;x 

i=i L4 J 

y[3](x;x) = y[2](x; x) + ]T [7(20*? - 3^)<y(x)(2(kf - 3^)>, lX/aiS 

i=l 

+ 15(12^ - l)^(y(x)(12^ 2 - 1)^)^x1255 

with i! = i(mod 2) + 1. Consequently, we have explicitly derived the exact expression of 
the macroscopic structures of fluctuating systems in arbitrary dimensions. Our result is 
generally applicable for any continuous systems regardless of the shape of the macroscopic 
structures. No more numerical fitting is needed. Besides, the obtained result can greatly 
help further analytical investigation of fluctuating systems. 



III. THE DETRENDED VARIANCES AND THE CORRELATION FUNCTION 



In the following, we would like to obtain the explicit and exact relations between the 
detrended variance functions and the correlation function. It will then facilitate the scaling 
analysis, which will be undertaken in the next section. First, the correlation function is 
defined as 

G(r) = [y(x)-y(x + r)] 2 (11) 

with the overbar denoting the average over the whole system. With some calculation, it is 
straightforward to obtain the relation between the correlation function G(r) and the average 



of the original variance function w 2 (J\f =1 If, x) (i.e., taking the average over all the windows 

of size rnU k)- 

/ d \ d T 9 rk 1 1 

(12) 
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G(r). 



Next, the average of the [gj-th degree detrended variance function is related to the average 
of the original variance function as: 
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Subsequently, 
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G(r). (14) 
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By employing the properties that the Legendre polynomials are either even or odd functions 
of their arguments and the correlation function G(r) is an even function of all r iy we obtain 
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G(r). (15) 



Substituting Eqs. (fT2"|) and (|T5|) into Eq. (|T3*j) . we consequently obtain the exact relation 
between the correlation function and the average of the [g]-th degree detrended variance 



function as follows 
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with the kernel given by 

K[q](r)= E (nf2(2n, + l) 

IWI<<? "=i L 



(17) 



As an illustration, let us consider the case with d = 2. We explicitly list out fT[0](r) to 
A"[3](r): 



K[0](t) = H(2l t -2r t ), 

i=i 

K[l](r) = K[0}(r) + £ - 6r, + 4?pJ (2l v - 2r t ,), 

K[2](r) = K[l](r) + £ - 10r, + 20^ - 12^ {21, - 2r v ) + J[ \2l t - 6r, + 4- 
tf[3](r) = tf[2](r)+£ 



2 / ' 
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+ [2k- 10n + 20^ - 12^ j Ul v - 6rv + 4^| 



with %' = i(mod 2) + 1. In a word, the correlation function and the [q]-th degree detrended 
variance function are closely related through a kernel function. Thus, before undertaking 
the scaling analysis of the fluctuating systems, one should first investigate the properties of 
these kernel functions. 

IV. SCALING ANALYSIS 



In this section, we will first explore some important properties of the kernel functions 



and then analyze the scaling relations among w 2 (Y\f =1 k; x), w 2 [q\(Y\f =1 li;Sc), and G(r 
The kernel K[q](r) under consideration can be recast as K[q](r) = J2\{ n }\< q Ui=iF ni (ri) 
By some calculation, we first obtain the following property of F n . (r-j): with being any 
nonnegative number (including non-integers), 
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We then will employ the above relation to compute Jj-jd ^ Q ; ^ K[q](r)r 2l3 dr (with r = |r|), 
in order to explore the scaling relations between the correlation function and the detrended 
variance functions. For (3 being a positive integer, 
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m, 

otherwise. 



1, 2, . . . , 7i j 



1; 



(2m.i+l)r(m l -ni+l)r(m i +n i +2) ' 

With some calculation and by induction, it can be shown that the above integral is equal 
to for (3 = 1, 2, . . . , q. In addition, in many physical situations, one usually takes the 
observation window with equal side lengths (i.e., = I for all i) and, then, the above 
integral is proportional to / 2 ^+ 2d ; for = q + l,q + 2, .... 

If j3 is not an integer, the multinomial expansion is no more applicable. The calculation 
becomes extremely complicated and no general closed form can be attained. Since most 
physical situations are involved with the substrate dimension d = 2, we then choose to 
undertake the calculation for d = 2 with the side lengths of the observation window li = 
1% — I. First, consider the integral 

H /3 (l,n 1 ,n 2 ) 



Ixl 



7711+712+2/3+2 



2(m + n 2 + 2(3 + 2) ^ 
with the incomplete Beta function 16] 
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(v,H;x) = [ X t^il-tY^dt. 
Jo 



By employing the above integral and the expansion formula for the Legendre polynomial, 
we then derive 



/ rfrF ni (r 1 )F n2 (r 2 )(r 1 2 + r 
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-4( 2ni + 1) t ^ML^ H^j + 1,0) (20) 

with (n)j = n(n — 1) • • • (n — j + 1). Subsequently, by using Eqs. (jTUjl and (|2(jp and the 
following relations for the incomplete Beta function |l6|| 

B(p, v, x) = T ^ T +"1 ~ B <". K 1 " x) (21) 



we finally obtain 

ir[g](r)(r? + r^)^dr 

J2 f dr x f dr 2 F ni {r x )F n2 {r 2 ){r{ + riy 
7~^„Jo Jo 



ni+n,2<q ' 



w r2v +2g+2 )g (-p) ff.tn. E (2ni+1)(2 „ 2+1) 



oo oo 

xEE 

j=0 fc=0 



(-ni)j(l + ni)j(-n 2 )fc(l + n 2 ) fc 
(2j + l)(2fc+l)(j + A; + /3 + 2)(jU!) 2 



1 

+ 



(A; + (j + I] 



2 



x yWf (2n | nfo n | n y (~n) fc (l + n) fc 
^ ^ i (2n + 1)(g n + 1) t (2fc + l)(2^ + 2/5 + 3)(^ 



_ 2 /3+ 



+ 



1 



.(A: + lUi (l/2)„ + iJJ 
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Next, we will proceed to analyze the scaling relations among the correlation function, 
the average of the original variance function, and the averages of the detrended variance 
functions. For a fluctuating system consisting of a macroscopic structure, we can recast 
y(x) as y(x) = y m acro(x) + y sto (x). Here, y m acro(x) represents the macroscopic structure, 
which is expected to be continuous and smooth. |/ s to( x ) represents the stochasticity relative 
to the macroscopic structure. From the definition of the correlation function G(r), G(r) can 
also be separated into two parts G m acro(r) + (jr sto (r). It is physically reasonable to assume 
that G m acro( r ) is analytic and thus can be expressed in terms of a power series expansion. 
Since the correlation function is a "difference" correlation function, we further expect that 
G , macro (r) is an even function of all n and thus G ma cro(r) = Ej{ n }|=i &{n} Ilf=i^ n< - From 
Eq. ()12|). it is straightforward to see that the average of the square of the original variance 
function will bear exactly the same scaling behavior as the correlation function. In most 



physical situations, the dominant term of G sto (r) is proportional to r 2a with a characteristic 
scaling exponent a, which is usually not an integer. Thus, for the experimental systems 
with the formation of the macroscopic structure, it is hard to measure the scaling exponent 
a directly, if r 2a is not the leading term in the correlation function. By substituting the 
above expansion of G(r) into Eq. and employing the properties of the kernel functions 
(which we just derived in this section), we rigorously show that, by raising the degree 
of the detrended variance function, the contribution due to macroscopic structure can be 
successfully suppressed. The above result is applicable in any dimensions. In addition, at 
least in two dimensions, we rigorously show that the detrended variance functions do retain 
the scaling exponent a (if it is not an integer) originated from the stochastic nature of the 
fluctuating systems. 



V. APPLICATION IN INTERFACE SUPERROUGHENING 

Dynamics of superrough growth processes can be modeled by stochastic partial differential 
equations: 

^^ = F({Vh}) + V ^t) + v (24) 

with h(x,t) denoting the interfacial heights, x being a position vector in d dimensional 
substrate, F describing the surface relaxation mechanism, rj representing the noise, and Vq 
denoting the average interface advancing velocity. The global fluctuation of the interface 
can be measured by the global interfacial width w(L,t) (with L denoting the linear size of 
the substrate), which is scaled as 

f t a l\ for t < L z - 
w(L,t)~{ (25) 

[ L a , for t > L z . 

Here, the two independent scaling exponents a and z are given the names the global rough- 
ness exponent and the dynamic exponent, respectively. The superrough interface is defined 
by the global roughness exponent a > 1, i.e., w sat urated(£)/£ — > oo as L — > oo. Qualitatively, 
the superrough growth processes gradually develop global mountains or valleys as time in- 
creases. By taking a snapshot of (d + l)-dimensional superrough interface morphology at a 
certain time, the formalism we just developed in the previous sections can be applied right 
away. Besides, one of the most distinct features of superrough growth processes is that the 
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local scaling differs from the global scaling; namely, the local interfacial width 



w(l,t) ~ < 



t a/z , for t < l z ; 

t (a-a loc )/ Zl a lo% for p < f ^ £*. ( 2 6) 
L a-a loc p lo for ^ > £ 2 



with / denoting the linear size of the local observation window. The exponent a\ oc is usually 
called the local roughness exponent. 

In the following, we will discuss two important examples in interfacial superroughening 
processes as applications. The first example is the experiment of nickel surfaces grown 
by pulse-current electrochemical deposition in 2+1 dimensions (J). The surface morphology 
taken by atomic force microscopy clearly shows the formation of large mounds. The values 
of the scaling exponents, measured in experiments Q], are a = 2.8, ot\ oc = 1.0, and z = 4.1. 
It is known that electrochemical noise can be regarded as temporally correlated noise [l7 1 . 
Thus, we propose the following stochastic partial differential equation to describe this growth 
process: 

r>h(v A 

= -z/V 4 /i + r](x,t) +v (27) 



Of 

with the noise correlation 



(r?(x, t)77(x\ t')) = D\t- t'^-^x - x'). (28) 

In the RHS of Eq. (|2T|). the first term accounts for the surface diffusion of deposits to 
the position of lowest chemical potential and the second term accounts for the temporally 
correlated noise. In 2+1 dimensions, by simple scaling anal ysis , we obtain a = 1 + 46 
and z — 4. In addition, generalizing the analysis in Ref. 18j to 2+1 dimensions, we 
obtain a loc = l[20j. Thus, the pulse-current electrochemical deposition processes can be 
well described by Eq. (12 7|) with the noise correlation index 9 = 0.45. 

Next, let us consider a broader class of growth processes in d + 1 dimensions: 

dh(x, t) = ( _ 1)m+Vv2m/l + ^ t) + Vq (2Q) 

with m = 1,2,3,.... The noise can be either the white noise (77 (x, t)rj(x.', £')) = DS(t — 
t')5(x — x') or the correlated noise (r]('x,t)r](yi , ,t')) = D\t — t'| 26,_1 |x — x'| 2p ~ d . The noise 
correlation indices p and 9 are restricted in the range (0, 1/2) to avoid the unphysical di- 
vergence. Generalizing the analysis in Refs. [l8j and 19( to d + 1 dimensions, we obtain 

11 



(a, z) = (^y^, 2m) for the cases with white noise and (a, z) = + p + 2m0, 2m) for 

the cases with correlated noise. Moreover, for the above class of growth processes with 
a > 1, the local roughness exponent ai oc is always equal to 1, independent of d, m, p, and 
20]. Is this just a coincidence? The further analytical calculation shows that for the above 
mentioned growth processes with a > 1, G macro (r, t) (the part of the correlation function 
attributed to the macroscopic structure formation) takes the form Y^Li A2 q t 2( - a ~ q ^ z r 2q and 
G s t (r,t) ~ r 2a for r z ^ t ^ L z 20]. This result confirms our scaling arguments in the 
previous section. Note that the leading term of G(r, t) is t 2 ( a ~ 1 )/ 2 r 2 for r z <^ t <^ L z . Recall 
that the square of the the local interfacial width bears the same scaling behaviors as the 
correlation function. Consequently, we have ai oc = 1, an intriguing property of superrough 
interfaces. To successfully extract the scaling exponent a from the experimentally measured 
data of interface morphology, the detrended variance function of degree [q] with the integer 
q > a — 1 needs to be employed. 



VI. CONCLUSION 



In conclusion, fluctuating systems with the formation of macroscopic structures are com- 
monly observed in Nature, for example, growth by molecular beam epitaxy, growth by elec- 
trochemical deposition, cultivated brain tumor growth, spontaneous imbibition of viscous 
fluids in porous media, the mosaic structure of DNA sequences, physiology signals, and the 
clustering features of seismic sequences. The "generalized detrended fluctuation analysis" 
method is designed to determine the correct scaling behaviors of fluctuations in the presence 
of possible trends (or macroscopic structures) without knowing their origin and shape. In 
this paper, we undertake an extensive study on this method in arbitrary dimensions. We not 
only derive the explicit expression of the macroscopic structures, but also obtain the exact 
relations between the correlation function and the detrended variance functions. Through 
the rigorous analysis of the kernel functions, we explicitly show that the detrended variance 
functions can successfully suppress the influence of the macroscopic structures on the scal- 
ing behaviors and, at the same time, keep the scaling behaviors of fluctuations unaltered. 
Finally, as applications, we discuss some important examples in interfacial superroughening 
and give explanations about an intriguing property of superrough interfaces aq oc = 1- All 
the obtained results are exact and applicable for continuous fluctuating systems in arbitrary 
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dimensions. 



Acknowledgments 

The work of N.-N. Pang is supported in part by the National Science Council of the 
Republic of China under Grant No. NSC-94-2112-M002-021. The work of W.-J. Tzeng is 
supported in part by the National Science Council of the Republic of China under Grant 
No. NSC-94-2112-M032-008. 



[1] J. Krug, Adv. Phys. 46, 139 (1997). 

[2] T. Halpin-Healy and Y.-C. Zhang, Phys. Rep. 254, 215 (1995). 

[3] S. Das Sarma, C. J. Lanczycki, R. Kotlyar, and S. V. Ghaisas, Phys. Rev. E 53, 359 (1996). 

[4] M. Saitou, Phys. Rev. B 66, 073415 (2002). 

[5] A. Bru et al., Phys. Rev. Lett. 81, 4008 (1998). 

[6] J. Soriano et al., Phys. Rev. Lett. 95 104501 (2005). 

[7] H.-N. Yang, Y.-P. Zhao, G.-C. Wang, and T.-M. Lu, Phys. Rev. Lett. 76, 3774 (1996). 

[8] J. H. Jeffries, J. K. Zuo, and M. M. Craig, Phys. Rev. Lett. 76, 4931 (1996). 

[9] C.-K. Peng et al., Phys. Rev. E 49, 1685 (1994). 
[10] S. V. Buldyrev et al., Phys. Rev. E 51, 5084 (1995). 
[11] K. Hu et al., Phys. Rev. E 64, 011114 (2001). 
[12] A. Bunde et al., Phys. Rev. Lett. 85, 3736 (2000). 
[13] S. Bahar et al., Europhys. Lett. 56, 454 (2001). 

[14] L. Telsca, V. Cuomo, V. Lapenna, and M. Macchiato, Geophys. Res. Lett. 28, 4323 (2001). 
[15] N.-N. Pang and W.-J. Tzeng, Phys. Rev. E 69, 031108 (2004); Phys. Rev. E 70, 036115 
(2004). 

[16] See, for example, An Atlas of Functions, eds. J. Spanier and K. B. Oldham (Springer- Verlag, 
Berlin, 1987); it Integrals and Series, eds. A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev 
(Gordon and Breach, New York, 1986). 

[17] B. Roseler and C. A. Schiller, Materials and Corrosion 52, 413 (2001). 

[18] N.-N. Pang and W.-J. Tzeng, Phys. Rev. E 70, 011105 (2004). 



13 



[19] N.-N. Pang and W.-J. Tzeng, Int. J. Mod. Phys. B 26, 3429 (2001). 
[20] The detailed mathematical derivation will be published elsewhere. 



14 



